Mean Free Path and Energy Fluctuations in Quantum Chaotic Billiards 
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The elastic mean free path of carriers in a recently introduced model of quantum chaotic billiards in 
two and three dimensions is calculated. The model incorporates surface roughness at a microscopic 
scale by randomly choosing the atomic levels at the surface sites between — W/2 and W/2. Surface 
roughness yields a mean free path I that decreases as L/W 2 as W increases, L being the linear system 
size. But this diminution ceases when the surface layer begins to decouple from the bulk for large 
enough values of W , leaving more or less unperturbed states on the bulk. Consequently, the mean 
free path shows a minimum of about L/2 for W of the order of the band width. Energy fluctuations 
reflect the behavior of the mean free path. At small energy scales, strong level correlations manifest 
themselves by small values of T?{E) that are close to Random Matrix Theory (RMT) in all cases. 
At larger energy scales, fluctuations are below the logarithmic behavior of RMT for I > L, and 
above RMT value when I < L. 



PACS number(s): 72.15.Lh, 05.45. +b 
I. INTRODUCTION 

Quantum systems in which the scattering of carriers 
by bulk impurities or defects is negligible are tradition- 
ally considered as ballistic. An interesting example of 
these systems are the so called billiards: a region of space 
in which all scattering occurs at the boundary. The 
ballistic character of these systems is usually taken for 
granted, even in|-|the case that their classical analogs be- 
have chaoticallyu. This is so despite of the fact that in 
wave mechanics the localized character of the scattering 
centers does not guarantee that their effects be also lo- 
calized. On the other hand, experimental data for the 
mean free path of microcavities do always correspond to 
the bulk material from which the cavity was produced. 
For instance, Chang et aln give a transport mean free 
path of the bulk material which is ten times larger than 
the cavity diameter. 

In this work, we characterize quantum billiards by the 
standard magnitude used in transport studies, i.e., by 
its mean free path. Surprisingly, the calculated elastic 
mean free path I is smaller than the linear size of the 
system for physically relevant values of surface rough- 
ness. Afterwards, the consequences of such small val- 
ues of I are studied on measurable magnitudes like the 
long-range energy correlations within the spectrum. Our 
calculations are done for a model of quantum chaotic bil- 
liard recently proposed by the authorstm. The model 
considers a system which is perfect all throughout the 
bulk, whereas its surface shows roughness at a micro- 
scopic scale. Bulk disorder is not considered since, as re- 
marked above, experimental mean free paths of the bulk 
systems are usually much larger than cavity diameters. 
An efficient implementation of this idea was developed 
in our previous workB: a tight-binding Hamiltonian with 



a single atomic level per lattice site in which the ener- 
gies of the atomic orbitals at the surface sites are cho- 
sen at random. We proved that, in the thermodynamic 
limit, this model, shows the accepted signs of quantum 
chaotic behavioru~D: an energy spectrum characterized 
by a Wigner-Dyson distribution of nearest-level spac- 
ings, and typical scars in the probability amplitude of 
wavef unct ions . 

Section II introduces the quantum billiard model. The 
elastic mean free path is calculated in two complemen- 
tary ways in Section III. Long-range energy fluctuations, 
which arc known Act-be a clear hallmark of diffusive or 
ballistic behaviorsa Ej, are analyzed in Section IV. The 
paper ends with a Section that collects our main results. 



II. QUANTUM BILLIARD MODEL 

Our quantum billiard model is defined by the following 
tight-binding Hamiltonian: 
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where the operator c, destroys an electron on site i, and 
tij is the hopping integral between sites i and j (the 
symbol (ij) denotes that the sum is restricted to near- 
est neighbor sites). We take tij = t = — 1 and consider 
square and cubic lattices for 2D and 3D systems, respec- 
tively. The energy u>i of a given atomic level at the sur- 
face (S) site i, is randomly chosen between —W/2 and 
W/2, whereas other sites have a constant energy equal 
to zero. Calculations have been carried out for squares 
of side up to L = 130 and for cubes of side up to L = 20, 
where L is given in units of the lattice constant. The 
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Schwarz algorithm for symmetric band matrices^ was 
used to compute the whole spectrum. 

Note .-that this model, in contrast with standard 
billiardsO, has two length scales: the system size L and 
the lattice constant a. Thus, even in the macroscopic 
limit (L/a — > oo) microscopic roughness remains, and is 
felt by quantum particles, i.e., particles characterized by 
a wavelength of the order of a. 



III. ELASTIC MEAN FREE PATH 

The ratio l/L, I being the mean free path, is the pa- 
rameter that controls whether a system is, diffusive or 
ballistic and so the type of fluctuationsaliJ. /The, stan- 
dard definition of the elastic mean free path ial3L3: 
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where r k is the relaxation time and i>k the velocity of 
a state with momentum k. We calculate the selfenergy 
induced by disorder in an ordered state of momentum k 
and energy e k from Dyson equation: 



E(e k ) = <G(e k )> -G^k) 



(3) 



where E is the selfenergy operator, Go the unperturbed 
Green function, and (G) the perturbed Green function 
averaged over disorder realizations. Then, in order to 
calculate the matrix element E(k, k; e k ), we only have to 
determine the inverse of the averaged perturbed Green 
function in the k-space basis using its site representation: 

(k|(G(e k )f *|k) =^V'k(r 4 )V'k(r,)(G(z,j;e k ))- 1 , (4) 

where "0k (r^) are the wave functions of the ordered clus- 
ter: 
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y^sin(k x Xi)sin(k y yi)sin(k z Zi) . (5b) 



It should be noted that non-diagonal elements in 
Eq.(4) are essential in this case, due to the strong break- 
ing of translational invariance introduced by surface dis- 
order. On the other hand, we have checked that a calcu- 
lation of the selfenergy due to bulk disorder without the 
incorporation, of non-diagonal elements gives very accu- 
rate resultsE3lI3. 

The qualitative behavior of X, can be obtained by 
second-order perturbation theoryta. Averaged real-space 
matrix elements of the selfenergy operator are approxi- 
mately given by: 
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(6) 



for those sites i for which the diagonal energy u>i fluc- 
tuates, i.e., for the surface sites. Using the basis given 
by Eq.(5) to transform to k-space and assuming a con- 
stant density of states (l/(4d), being d the dimension 
of the physical space) to get a rough approximation for 
the imaginary part of the diagonal elements of the un- 
perturbed Green function (Im[Go(i, i; e k )] ~ 7r/(4d)), we 
obtain: 
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Substituting the selfenergy in Eq.(g) for this approximate 
value and using the 2D dispersion relation of the band 
e k = — 2[cos(fc :c ) + cos(fcy)]) to get an explicit equation 
for fi.v k , we finally have: 
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The final expression for I only uses dimensionless mag- 
nitudes, i.e., disorder energies W are measured in units 
of the hopping energy t, lengths in units of the lattice 
constant a, and k in units of a -1 . The corresponding ex- 
pression for 3D has an additional term within the square 
root. 

When second-order perturbation theory is used to get 
an approximate value of the mean free path induced by 
bulk disorder, a rather good approximationt5ll2l is ob- 
tained due to the homogeneous character of both the 
perturbation and the real-space Green function. Nev- 
ertheless, when the same scheme is applied to billiards, 
it gives the correct L and W~ 2 dependences but fails 
to describe the variation of I within the band (see our 
numerical results shown in Figj^). 

Comparison of the mean free path of a billiard (Eq. (@)) 
with the result obtained for bulk disorderE3 reveals an 
essential feature of the present model: whereas for bulk 
disorder I is independent of the system size L, in the 
present case it depends linearly on L. This result is a 
consequence of having the scattering centers located at 
the surface, an essential feature of billiards. Therefore, 
ballistic (/ > L) or even clean behavior (I ^> L) could 
naively be expected over the whole band for any value 
of W. Our results will show that this is not the case for 
physically relevant values of W. 

Fig.|l| shows typical numerical results for the mean free 
path in 2D and 3D. Note that I is indeed proportional 
to the system size L in agreement with the approximate 
result reported in Eq.(||). The dependence of the mean 
free path on the energy reflects the nature of the eigen- 
values discussed in our previous worksau. In 2D, I in- 
creases near the edges and at the center of the band, 
due to ±he presence of quasi-ideal states in these energy 
regionsM. Quasi-ideal states do also appear in three di- 
mensions near the band edges and this fact is also re- 
flected in the behavior of the mean free path. The energy 
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dependence of I is significantly weaker than in the case of 
bulk disorder. Actually, I is almost constant except near 
the band edges and the band center. 

Eq.(||) shows that the mean free path monotonically 
decreases with W. However, we note that the definition 
of the mean free path (Eq.(g)) is intrinsically perturba- 
tive. As a consequence, the results for large W cannot 
be trusted. In fact in the limit of large W the surface 
decouples from the bulk, leaving a fully ordered cluster 
of linear size L — 2. Consequently, the mean free path of 
the inner part should tend to infinity. In order to handle 
this problem, we have calculated the mean free path in 
the large W limit by considering the perturbation of the 
ordered states of the (L — 2) d cluster due to the addition 
of the surface layer. 

Firstly, the averaged Green function of the whole clus- 
ter of lineal size L is calculated for the eigenenergies of 
the smaller cluster of linear size L — 2. Secondly, the 
averaged Green function of the inner cluster is obtained 
using a projection operator P, and finally, the selfenergy 
operator is defined by means of an equation quite similar 
toEq.(|): 



E(e k ) = [p+(G(e k )}p| - G^(e k 



(9) 



Since Eq.(^|) is an operator definition, any basis can 
be used in the numerical calculation. We use the tight- 
binding basis in which the projector P is particularly 
simple. Thus, the diagonal elements of the selfenergy in 
the k-space basis are given by: 



(k|S(e k )|k) =^^. k (r i )V; k (r i )S(i,i;e k ) 



(10) 



of long-range correlations within the energy spectrum. 
Our data show that the scaling behavior given in Fig.[| for 
two dimensional billiards holds for billiards in spaces of 
larger dimension. This result is not surprising according 
to well-known results for bulk disorder: mean free path 
is weakly dimension .dependent in those regions in which 
it is almost constantEj. 

Let us mention that although nearest-neighbors spac- 
ing statistics has been proved to follow Wigner-Dyson 
statistics in our previous worktj, finite clusters only show 
"nice" GOE statistics for a range of W values close to 
the bandwidth. Fig.|| illustrates this fact. It shows the 
variance of nearest-neighbor spacing distribution as a 
function of W for several system sizes. As expected, 
chaotic behavior is not found when disorder is either very 
strong or very weak. Beyond a value of W ~ 10 which 
is of the order of the band width, the variance becomes 
larger than the value corresponding to the Wigner-Dyson 
statistics (0.286), indicating that the surface starts to de- 
couple from the bulk. For this value of W, most of the 
energy levels which are mainly localized on surface sites 
lie outside the band and follow Poisson statistics. Also 
for small W values, the value becomes larger than 0.286. 
This indicates a clean behavior of our system: statistics 
of clean (orderedjusystems is more or less Poisson up to 
some energy scalaLa. The dependence of the variance on 
the system size is also illustrated in Fig.||. When L in- 
creases from 32 to 90, the lower value of W below which 
the variance is no longer that of GOE decreases. Fur- 
ther increases of L would reduce to zero the widths of 
the energy regions where quasi-ideal states appear, and 
thus, the variance would become that of GOE through 
the whole band, no matter how small W would be. 



This alternative perturbation theory provides us with 
a second value of the mean free path. Of course, the 
larger value is the only one which is physically relevant 
(perturbation theory logically starts from the less per- 
turbed system) . Although we have not found an analytic 
approximation for this second value of I, our numerical 
results show that I is indeed linearly dependent on L (see 
Fig.^). Apart from some features that can be related to 
singularities of the density of states, this new value of I 
behaves in the same way as the standard value: it shows 
considerable fluctuations around a mean value that is al- 
most constant within the band. Similar results have been 
found for three dimensional billiards. 

Fig.|| shows the scaling behavior of both values of I 
with disorder. The relevant (larger) values of I closely 
follow an W~ 2 law for small W values while the scaling 
is close to W for large W values. These scaling laws do 
not depend on the system size. On the contrary, adding 
the external shell to the cluster is a perturbation that 
decreases in importance as the cluster increases (see the 
behavior of triangles in the left part of Fig||) . Our results 
show that the ratio l/L can reach values as small as 0.4 
close to W = 6. We will see in the next section, that 
the value of this ratio strongly determines the behavior 



IV. ENERGY FLUCTUATIONS 

In this section we analyze the relationship between 
mean free path and measurable physical magnitudes. 
Although conductance was our first candidate, finally 
we chose long-range energy fluctuations due to the im- 
portance they have in analyzing quantum chaotic sys- 
tems. Moreover, while the behavior of energy fluctua- 
tions within the spectrum is an intrinsic property, con- 
ductance depends on the way the system is attached to 
the metallic leads. 

The variance of the number of levels in a given en- 
ergy interval is a well-known statistics usually referred 
as Y?{E) statistics!! The standard way of measuring 
E 2 (E) proceeds firstly to transform real spectra onto un- 
folded spectra: each real spectrum {e^} is mapped onto 
an unfolded spectrum {Ei} through Ei = N(ei) ) where 
N(e) is the averaged number of levels up to an energy e. 
Secondly, the variance of the number of levels found in 
an interval of fixed length E is directly obtained from the 
unfolded spectra. In this Section, we present results for 
the variance of the number of levels obtained by a dif- 
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ferent numerical procedure. Of course, we have checked 
that overall features are procedure independent, i.e., our 
method coincides with the standard one when subtleties 
are ignored. Nevertheless, fine details of the number vari- 
ance are better given by a more direct definition. 

Let us sketch the numerical procedure. Some energy 
interval [ei, 62] of the band spectrum is chosen. The mean 
number of levels within the interval is given by: 



N=(N(e 2 )-N(e 1 ) 



(11) 



where N(e) gives the total number of states below energy 
e for a generic disorder realization and averaging over 
disorder configurations is done. In the same way, the 
mean of the squared number of levels is given by: 



N* = / [JV( ea ) - iV(ei)] 1 



(12) 



The variance of the number of levels contained in the 
energy interval [£1,62] is simply: 



£ 2 (A) = N 2 — N 



(13) 



In this way, just a value of the T, 2 (E) function is obtained. 
This sequence is repeated a large number of times for ran- 
domly selected energy intervals until a relatively smooth 
value is obtained for the variance of the number of levels 
averaged over energy intervals containing the same num- 
ber of levels. The last step implies averaging the number 
variance over some selected region of the energy band. 

Note that only fluctuations induced by disorder are 
taken into account by our method. Therefore, number 
fluctuation is strictly zero for an isolated spectrum. On 
the other hand, the standard method can be used even 
in this situation. The variance is defined by: 



Y, 2 (E) = [N(Eo + E) - N(Eo) - Ef 



(14) 



The subtlety within this definition is the way followed 
for unfolding just one spectrum. Sometimes, the asym- 
totic form of iV 6 _ >00 (e) (Weyl formula)!! is used even in 
the lower part of the spectrum and number fluctuations 
are measured relative to N e - tOQ (e). Most frequently, a 
large number of levels is used to measure the average 
level spacing and then fluctuations relative to the mean 
value are measured in smaller intervals. In this way, it 
has been proved that crystalline spectra are uncorrelated 
to some energy extent, and consequently, Y?(E) statistics 
is Poissonlij: 



^Poisson (E) 



E 



(15) 



Let us remark that our procedure is quite appropriate 
for billiard ensembles as it is our case (surface roughness 
is random and several samples can be generated for a 
given value of surface disorder W). On the other hand, 
it gives no fluctuations for conventional billiards (just one 
spectrum) . 



We have checked that our numerical procedure recov- 
ers both the analytical results for random series of energy 
levels (see Eq.(|l^)) and for the eigenvalue series of ma- 
trices belonging to the Gaussian Orthogonal Ensemble 
(GOE)&3: 



J GOE 



(25) « -5 



ln(27r£;) + 1.5772 



(16) 



The logarithmic dependence on E makes this value much 
smaller than the one corresponding to uncorrelated spec- 
tra. The well-known spectral rigidity of GOE follows: 
the number of levels in any energy interval is very close 
to its mean value. 

As said in the last part of the previous Section, we 
routinely calculate the variance of the nearest-level spac- 
ing distribution in order to check that the Wigner-Dyson 
value of 0.286 is recovered by our numerical code. This is 
always the case for "reasonable" values of W (see Fig.^) 
in spite of dealing with systems of finite size. Notice 
that the calculation of S 2 (E) for very small values of W 
is difficulted by finite size oscillations in the density of 
statesO. 

Fig.pl shows the typical behavior of the number vari- 
ance Tr{E) in different transport regimes. Results cor- 
respond to the spectra of a large collection of 32 x 32 
clusters (one thousand samples). All states within the 
central part of the band (— 3 < e < 3) have been consid- 
ered. The black thin line corresponding to W = 2 shows 
a typical ballistic behavior: fluctuations are really small 
in a large energy interval (the number of levels fluctu- 
ates in less than one level for intervals containing 100 
levels). On the other hand, the black thick line shows a 
typical diffusive result that corresponds to W = 10. In 
this case, number fluctuation follows GOE statistics for a 
small number of energy levels (w 8), remains below GOE 
up to about 35 levels, and is well above GOE values at 
larger energy scales. In any case, fluctuations are still 
much smaller than Poisson values (£p oisson (100) = 100). 
The gray thick line shows results corresponding to an 
intermediate situation (W = 6). A natural physical ex- 
planation of such different behaviors is provided by the 
ratio of the mean free path to the cluster side: l/L is 
about 0.7 for the diffusive system but 2.6 for the ballistic 
one. Although both values are not too far from the ac- 
cepted assumption I ~ L, long-range energy correlations 
within the spectrum are quite different. Our results for 
the intermediate value of W prove that the transition 
from one regime to the other is smooth. 

Let us further comment our numerical result for bidi- 
mensional billiards of side L = 32 (FigJ|). Firstly, a 
continuous increase of the number of levels following 
GOE statistics is observed for increasing values of W 
(W = 2,6,10). Secondly, it follows a region in which 
Y?(E) deviates from GOE showing a saturation effect. 
And finally, an almost linear increase of the variance fol- 
lows at larger energies. Our qualitative interpretation 
of this behavior is the following. Strict GOE statistics is 
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limited by the maximum degeneracy of crystalline bands. 
For a 2D cluster of the square lattice with open boundary 
conditions, maximum degeneracy is 2 although splitting 
of the 8-fold levels corresponding to periodic boundary 
conditions is really small. Surface roughness produces a 
mixing of this degenerate states that is well described by 
the random matrix elements of matrices of the GOE. This 
suffices to give Wigner-Dyson statistics for the nearest- 
neighbor spacing distribution. When larger energy scales 
are considered, the fact that crystalline levels range from 
energy —4 to 4 implies an incomplete mixing of levels 
(in the limit of very small non-diagonal random cou- 
pling, T, 2 (E) would be constant beyond the number of 
levels corresponding to maximum degeneracy). This ex- 
plains saturation with variance values below GOE statis- 
tics. The almost linear Poisson-like increase of the num- 
ber variance at larger energy is a subtle effect associated 
to the randomness of the energy levels at the surface. 
Energy levels suffer uncorrelated shifts due to the forma- 
tion of surface resonances at uncorrelated energies. It is 
clear that for a large enough value of W, about 4L sur- 
face states appear outside the band. For the same rea- 
son, about AL surface resonances exist within the band 
for smaller W values. While saturation can be easily 
checked using matrices of random non-diagonal elements 
but fixed diagonal elements, the Poisson-like linear in- 
crease found at larger energies is a characteristic feature 
of our model. 

Similar results are obtained in 3D as shown by Figj^. 
A large collection (one hundred samples) of clusters with 
random surface roughness has been fully diagonalized 
and the number variance has been calculated. Fluctu- 
ations have been obtained for the central part of the 
band excluding the band center: —4 < e < —2. In this 
case, ballistic regime is clearly obtained for W — 2 (black 
thin line) whereas diffusive behavior is got for W = 10 
(black thick line). An intermediate behavior is obtained 
for W = 4 (gray thick line). Results can be interpreted 
as in 2D. Strict GOE statistics is always followed by at 
least six levels (typical 3D crsytalline degeneracy for open 
boundary conditions). A region showing saturation of 
S 2 (i?) follows (this is clearer for W = 2 than in the 
other two cases) . At larger energy scales an uncorrelated 
behavior of the variance can be seen. The larger slope 
of the linear region found for W = 10 as compared to 
W — 2, reflects the fact that surface resonances are bet- 
ter defined in the first case. About 125 energy levels 
closely follow GOE statistics for the intermediate degree 
of surface roughness (W = 4) . We believe that this prop- 
erty is related to a minimum value of the l/L ratio. The 
smallest value of the mean free path at fixed cluster size 
is obtained when the mixing of crstalline states is maxi- 
mum. Therefore, a larger number of energy levels shows 
correlations that are well described by RMT. Generally 
speaking, ballistic behavior (I > L) corresponds to a min- 
imum mixing of energy levels (only degenerate states are 
mixed by the surface random perturbation) while diffu- 
sive behavior (I < L) is seen when surface resonances 



dominate the value of the number variance. 

Fig.0 shows the increase in the number of strongly 
correlated energy levels for a fixed surface roughness 
(W — 2) as the size of the system increases. In this Fig- 
ure, a clear discontinuous change of slope can be seen for 
both curves. About 140 and 300 levels are strongly cor- 
related for L — 14 and 20, respectively, before the effect 
of uncorrelated surface resonances dominates. One could 
naively suggest an L 2 scaling for the number of corre- 
lated levels. Nevertheless, a precise assesment about the 
scaling of such critical energy would require further BHp 
merical simulations well beyond our actual possibilities!!^!. 

V. CONCLUDING REMARKS 

We have found that the mean free path in the model 
of quantum chaotic billiards discussed in this work can 
be substantially smaller than the system size. Actually, 
the ratio l/L reaches a minimum value of about 0.5 for 
the most relevant disorder values (W of the order of the 
band width). The transition from ballistic to diffusive 
behavior is governed by the disorder parameter W which 
simulates surface roughness at a microscopic scale. The 
behavior of the number variance statistics (Y?{E) statis- 
tics) clearly distinguishes between both regimes. After 
closely following RMT for a small number of energy lev- 
els, Y?{E) shows a tendency to saturation before increas- 
ing in an almost linear way at larger energies. The region 
of saturation extends up to a really large number of levels 
in the ballistic regime, while the linear increase quickly 
dominates in the diffusive regime. Therefore, Y?(E) is 
roughly below or above the values corresponding to GOE 
statistics for ballistic or diffusive regimes, respectively. In 
any case, both transport regimes are still characterized 
by having a mean free path of the order of the linear 
size of the system and Y?(E) values much smaller than 
E, i.e., much smaller than values corresponding to com- 
pletely uncorrelated energy spectra. The behavior of the 
number variance has been related to the structure of the 
particular type of random matrices produced by surface 
disorder. Preliminary studies show that these results are 
not restricted to the particular model analyzed in this 
work since similar results have been found for a model 
with real surface roughnesstH 
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FIG. 1. Mean free path (divided by the system size L) 
as a function of the energy (scaled by half the band width, 
4 and 6 in 2D and 3D, respectively) for the chaotic billiards 
discussed in this work. In the 2D case the results correspond 
to L—8 and 32 (circles and triangles) and W = 2, and in the 
3D case to L = 6 and 10 (circles and triangles) and W = 4. 
The fact that the data for the two sizes collapse indicates that 
I oc L (see text). 
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FIG. 2. Mean free path (divided by the system size L) as a 
function of the energy (scaled by half the band width) for the 
chaotic billiards discussed in this work. Results correspond 
to L—16 and 32 (circles and triangles) and W = 10. The fact 
that the data for the two sizes collapse indicates that I oc L 
(see text). 



FIG. 4. Variance of the nearest-level spacing distributions 
as a function of the disorder parameter W for the two (circles, 
L=32, triangles, L = 90) and three (squares, L=10) dimen- 
sional billiards here investigated. The results correspond to 
an average over an energy window of 0.2 around the energy 
— 1. The dashed horizontal line corresponds to the GOE vari- 
ance (0.286). 
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FIG. 3. Mean free path I divided by L as a function 
of W for the quantum chaotic billiards investigated in this 
work. Standard approach results are given by circles whereas 
our alternative calculation is represented by down triangles. 
Small symbols correspond to L = 16 and large symbols to 
L = 32. The results correspond to an average over an energy 
window of 0.2 around an energy equal to —2. 
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FIG. 5. T, 2 (E) as a function of E for the 2D billiard. The 
results correspond to L — 32 and disorder W — 2 (black thin 
line), W = 6 (gray thick line), and W = 10 (black thick line). 
One thousand samples and all states in the region — 3 < e < 3 
are included in the calculation. Circles show the variance 
value corresponding to GOE fluctuations. 
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FIG. 6. E 2 (£) as a function of E for the 3D billiard. The 
results correspond to L — 20 and disorder W = 2 (black thin 
line), W — 4 (gray thick line), and W = 10 (black thick line). 
One hundred samples and all states in the region —4 < e < — 2 
are included in the calculation. The dashed line shows the 
variance value corresponding to GOE fluctuations. 
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FIG. 7. Number variance as a function of the mean num- 
ber of levels of a 3D billiard of surface roughness given by 
W — 2 and two different sizes. Results for L = 14 are given by 
filled circles whereas results for L — 20 are given by empty cir- 
cles. The corresponding GOE values are shown by the dashed 
line. 
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